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Abstract 



The propagation of an initially localized excitation in ID incommensurate, 
quasiperiodic and random systems is investigated numerically. It is discovered 
that the time evolution of variances o~ 2 (t) of atom displacements depends 
on the initial condition. For the initial condition with nonzero momentum, 
o~ 2 (t) goes as t a with a = 1 and for incommensurate Frenkel-Kontorova 
(FK) model at V below and above V c respectively; and a = 1 for uniform, 
quasiperiodic and random chains. It is also found that a = 1 — (3 with 
(3 the exponent of distribution function of frequency at zero frequency, i.e., 
p(uS) ~ (as u — ► 0). For the initial condition with zero momentum, a = 
for all systems studied. The underlying physical meaning of this diffusive 
behavior is discussed. 
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I. INTRODUCTION 



The fast development of nanotechnology makes one- dimensional (ID) systems like quan- 
tum wireB and nanotubesll available in laboratories nowadays. The study of physical prop- 
erties such as localization, thermal conductivity, and electron conductance, etc., in these 
systems become more and more important!. In the past years, many works have been done 
on the waves (classical and quantum) propagation and localization in ID disordered media. 
Particular interests are paid to propagation of the classical waves in random media (see 
e.g.!), and electron transport in ID disordered solids!'!. In comparison with the disordered 
systems, much less is studied about the transport and diffusion properties of incommensu- 
rate and quasiperiodic systems, even though the incommensurate structures appear in many 
physical systems such as quasicrystals, two-dimensional electron systems, magnetic super- 
lattices, charge- density waves, organic conductors, and various atomic monolayers absorbed 
on crystalline substratesi. 

In this paper, we would like to study the classical transport of an initially localized 
excitation in ID incommensurate, quasiperiodic and random systems so as to better under- 
stand transport processes and relaxations properties in these systems. Recently, quantum 
diffusion in two families of ID models have attracted much attentio Typical examples 

are the kicked rotator and kicked Harper models from the field of quantum chaos, and the 
Harper model and the tight-binding model associated with quasiperiodic sequences. Many 
interesting dynamical behaviors, such as quantum localization and anomalous diffusion, and 
their relationships with energy spectra have been investigated in these systems! - !^. How- 
ever, the classical transport in incommensurate and quasiperiodic sytems and its relation 
with phonon frequency distribution, as well as the diffusive behavior dependence on initial 
condition, etc. have not yet fully investigated up to now. 

The information of disturbance spreading in such system reflects the interior structures 
of the underlying system. As we shall see later, the spreading properties are determined 
largely by the density of states, in particular by the phonon model near the zero frequency. 



II. MODELS AND NUMERICAL RESULTS 

A. Incommensurate chain 

The Frenkel-Kontorova (FK) mo del0 is invoked as a prototype of an incommensurate 
chain in this paper. This model is a ID atom chain with an elastic nearest neighbor interac- 
tion and subjected to an external periodic potential. Most works on this model in the past 
two decades have been concentrated on ground state properties and phonon spectra, etc. 
The ID FK model is described by a dimensionless Hamiltonian 
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were p n and x n are momentum and position of the nth atom, respectively. V is the coupling 
constant, and a is the distance between consecutive atoms without external potential. Aubry 
and Le Daeron0 showed that the ground state configuration is commensurate when a/27r 
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is rational and incommensurate when a/2n is irrational. For an incommensurate ground 
state, there are two different configurations separated by the so-called transition by breaking 
of analyticity predicted by Aubryll^. This transition survives the quantum fluctuationEJ. 
Moreover, in contrast to other ID nonintegrable systems such as the Fermi-Pasta-Ulam 
chain0, the FK chain shows a normal thermal conductivity^. For each irrational number 
a/2n there exists a critical value V c separating the sliding state (V < V c ) from pinned state 
(V > V c ). The V c = 0.9716354 ■ • • corresponds to the most irrational number, golden mean 
value a/2n = (y/E— l)/2. Without loss of generality, we restrict ourselves to this particular 
value of a in the numerical calculations throughout the paper, and it is approximated by a 
converging series of truncated fraction: F n /F n+ i, where {F n } is the Fibonacci sequence. 
The equation of motion for the nth atom in the FK model around its equilibrium position 

is 

^ = t/W! + Vn-i " [2 + V cos(x°)]^ n , (2) 

where x° n is the equilibrium position of the nth atom at ground state, and ip n is the normalized 
displacement from the equilibrium position. In fact, to obtain Eq.(H), we have written the 
whole displacement of particle as x n = x° n + eip n , where e(^C 1) is a small parameter. To 
quantify the disturbance spreading, the variance of displacements 

1 N 

^)=vEW*)-^(°)l J ( 3 ) 

iV n=l 

is calculated by two numerical methods. The first one is the Runge-Kutta method of the 
fourth order to integrate Eq. (^) for a given initial condition with free boundary. The second 
one is to find eigenfrequencies ujj and eigenvectors a n (j) of equation 

- uo 2 ^ n = V>n+i + V>n-i -[2 + V cos(a£)]V> n . (4) 
The solution of Eq. can be expressed in the following form: 

N 

= i A i cos{uJjt) + Bj sxa(ujjt)) a n (j) (5) 
i=i 

where the coefficients Aj and Bj are determined by initial conditions. Contrasting to the 
quantum diffusion, the classical evolution Eq. (|2|) is of the second order of derivative. Thus 
initial conditions for ip n and dip n /dt are needed. One of our main findings is that the 
spreading behavior depends on the initial condition. For the initial condition 



ip n = and dip n /dt = 8, 



n, no 



which is called type I, we have a 2 ~ t a , and a is equal to 1 and for V < V c and V > V c , 
respectively. For the initial condition 

ipn = <Wo and dipjdt = 

which is called type II, a 2 ~ t° for any V. (Of course, there is another type of initial 
condition, i.e., ip n = 5 ni „ and dip n /dt = 5 njTV0 . Our numerical calculations show that the 
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spreading behavior in this case is the same as that of type-I initial condition.) Figure [I] 
shows the typical time evolution <J 2 (t) for the FK model. In numerical calculations, we 
first obtain the ground state positions of N atoms in the FK chain by the gradient method 
for free boundary i.e., x = and = Na. The results of Fig. [I] are obtained by the 
integration method for N = F±g = 10946. <J 2 (t) is also calculated by the second numerical 
method for the FK chains of small size, which gives rise to the same results. 

It is worth pointing out that the above-mentioned results valid only for the evolution 
time less than a critical value t*, where t* ~ N/2v and v is the velocity of sound. For our 
FK model, v ~ 1. After this critical time, i.e., t > f, the power relation of <J 2 (t) ~ t a is 
destroyed due to the finite size effect. 

To get a clear picture of the spreading of a disturbance in an incommensurate chain with 
two different initial conditions, we plot ip n {t) in Fig. 2 for the FK chain of V = 0.4. The 
intensity of gray scale represents the amplitude of the displacement of the particle. Because 
of the huge amount of the data, we record ip n (t) at a time interval of 20 time steps, which 
leads to some discontinuity. Figure 2(a) demonstrates the evolution with initial condition 
ip n — and dijjn/dt = 5„ 5 „ , and Fig. 2(b) shows that with initial condition ip n = 5 ntno and 
dip n /dt = 0. The difference is clear. In the later case the disturbance spreads out in both 
direction, and the particle remains almost at rest after the disturbance passes it. However, 
in the first case, wherever the disturbance spreads, the particle will be excited and keeps 
moving. In the cantorus regime (V > V c ), the disturbance spreading in the FK chain is 
similar to the case in Fig. 2(b) regardless of the initial condition. 



B. Quasiperiodic and random chains 

We turn now to study of disturbance spreading in uniform, quasiperiodic, and random 
chains. The equation of motion can be written as 

- = k n ip n+1 + fc n _i?/> n _i - (k n + k n -i)%l) n . (6) 

If k n = k for all n, it corresponds to a uniform chain. For quasiperiodic chains, k n takes 
two values fci and k<i which are arranged according to some deterministic quasiperiodic sub- 
stitute rulesiZl. Here we discuss four types of quasiperiodic chains. They are Fibonacci, 
Thue-Morse, Rudin-Shapiro, and period-doubling chains, respectively. The substitute rules 
for them are: k\ — * k\kiik,i —>■ k\ (Fibonacci); k\ — > kik2,k2 — > kvky (Thue-Morse); 
k\k\ — ► k\k\k\k2-,k\k2 — > kik^ki, k 2 k\ — ► fc 2 fc 2 fci&2, k 2 k 2 — > A^A^A^i (Rudin-Shapiro); 
k\ — > kik2,k2 — > k%ki (period-doubling). According to the classification based on the eigen- 
values of generating matrix defined by Luck@, they are bounded (Fibonacci, Thue-Morse), 
unbounded (Rudin-Shapiro), and marginal (period-doubling). For comparison, u 2 {t) for 
random chain is also studied. In this case, the values of k n are taken k\ and k-i with the 
same probability. Figure |^ shows typical time evolutions of variances for four quasiperiodic 
and random chains. The disturbance spreading behaviors in these chains are the same as 
that of the incommensurate FK model at V < V c , namely, <J 2 (t) ~ t for all these chains with 
the initial condition of nonzero momentum [dip n (0)/dt ^ 0], and t° for all these chains with 
initial condition of zero momentum [dip n (0)/dt = 0]. 
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C. Relationship with phonon spectrum 

Figures 1-3 are the main results. They demonstrate that the disturbance spreading 
depends crucially on the initial condition. In the following, we would like to understand this 
peculiar behavior in terms of the phonon spectra. 

The coefficients Aj and Bj in Eq. (5) are: Aj = and Bj = a no (j)/uj for type-I 
boundary condition; Aj = a no (j) and Bj = for type II. Therefore the solutions of Eq. (0) 
are 

N 

3=1 
N 

ij) n = ^cos{ujjt)a n {j)a no {j) : type II, (7) 

3=1 

respectively. As N — > oo, we have 

T7 £ ~ / sin 2 (wt)a 2 o pM^/u; 2 , type I, 

1 9 f^max 

T7 X, V'n ~ / cos ( ut ) a n P(^)d^ type n > ( 8 ) 

respectively, where uj max /uj m i n is the maximum/minimal frequency of phonon spectra and 
p(u;) is the density of the phonon spectra. 

The difference between the integrands in Eq. (||) for type I and type II lies in the factor 
1/u 2 . As time increases, the dominant contribution of the integral in Eq. (§) for type I 
comes from the integrand around u = 0. The integrand for type II is an oscillated function 
of time t, and so is the integral. Therefore the reason for different behaviors of these chains 
for different initial conditions is due to the coefficients Bj in Eq. (|5|). If Bj is equal to zero, 
i.e., the initial condition with zero momentum, u 2 {t) is an oscillated function of time. If Bj 
is nonzero, i.e., the initial condition with nonzero momentum, <J 2 (t) is proportional to t. 

In fact, the integral in Eq. (|8]) for type-I boundary condition can be written as 

/ ^ tsin 2 (^)a 2 p(-)^/^ 2 . (9) 

If the distribution function of frequency has the scaling behavior p{ui) ~ u 13 at low frequency 
(u — > 0), then one has a = 1 — (3. For the uniform chain, it is well known that p{uo) = 
2/ {-k^uj^ — uj 2 ), thus (3 = 0. We also discover that for the Fibonacci chain and random 
chain, the distribution function of frequency at low frequency are the same as that of uniform 
chaini'El. To demonstrate this, we calculate the integrated distribution function of frequency 
(IDFF) for these quasiperiodic chains by directly diagonalizing chains of finite length, and 
plot them in Fig. |as a function of uj. The results suggest that for all these quasiperiodic 
chains the IDFF is proportional to u at low frequency, thus (3 = 0, which is the same as that 
of uniform and random chains so that the relationship a = 1 — (3 is satisfied for all these 
systems. 
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For the incommensurate FK model, it is well known that there is a zero- frequency phonon 
mode for V < V c , whereas there is a phonon gap for V > V c . From above discussion, we 
know that the zero-frequency phonon mode plays a key role in the time behavior of u 2 {t). 
The time behavior of <J 2 (t) in the incommensurate FK chain at V < V c suggests that the 
low frequency behavior of distribution function of frequency is the same as that of those 
chains discussed above. The curves shown in Fig. f| indeed demonstrate this. But for the 
incommensurate FK chain at V > V c , 0J m i n > 0, thus the integral in Eq. (|8]) for type I is an 
oscillated function and the time behavior of <J 2 (t) is also an oscillated function of time. The 
case V — V c is critical. The phonon spectrum of the FK chain at V c is different from that of 

V < V c . It has self-similar structure and is point spectrum [see Fig. 0(a)] . Therefore there 
is no inverse power relation between p(u) and u at low frequency. It implies that the results 
would depend on the length of the chain in numerical calculation. This is illustrated in Fig. 
|5|(b), where we plot u 2 {t) as a function of t for the FK chains of different length at V = V c . 

III. CONCLUSION AND DISCUSSIONS 

We have studied the disturbance spreading in incommensurate, uniform, quasiperiodic 
and random chains. We have found that the time evolution of variance u 2 {t) depends on 
the initial conditions. Its behavior is determined by the density of phonon frequency around 
zero frequency. For the initial condition of zero momentum, u 2 {t) ~ t° for all kinds of chains 
studied in this letter. For the initial condition of nonzero momentum, u 2 {t) ~ t a , a = 1 for 
uniform, quasiperiodic, random chains, and incommensurate FK chain at V < V c . Although 
other physical properties differs from system to system, the time behavior of <J 2 (t) are the 
same for all these systems. For the incommensurate FK chain at V > V c , cr 2 (t) is an oscil- 
lated function of time. This different behavior of the incommensurate FK chain at different 

V regimes might provide us a different approach to detect the transition by breaking of 
analyticity experimentally. 

This work was supported in part by grants from the Hong Kong Research Grants Council 
(RGC) and the Hong Kong Baptist University Faculty Research Grant (FRG). Tong's work 
was also supported by Natural Science Foundation of Jiangsu Province and Natural Science 
Foundation of Jiangsu Education Committee, PR China. 
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FIGURES 



FIG. 1. a 2 (t) of the incommensurate FK chains with different values of V for different initial 
conditions. The solid and dotted lines correspond to V = 0.4 and V = 1.6, respectively, (a) Initial 
condition ipn(0) = and dip n (0) /dt = <5 njno ; (b) initial condition ip = S n ^ no and dtp n (0)/dt = 0. The 
length of the FK chain is TV = 10946. Initial excitation locates at the center of the chain, namely, 
n = N/2 + 1. 

FIG. 2. Time evolution of the displacement ijj n (t) for an incommensurate FK chain of V = 0.4 
with different initial conditions: (a) Y> n (0) = and dip n (0)/dt = <5„ jno ; (b) tp = 8 njno and 
dip n (0)/dt = 0. The length of the FK chain is N = 1597. The initial excitation locates at 
the middle of the chain, i.e., no = 799. 

FIG. 3. (J 2 (t) of quasiperiodic and random chains with different initial conditions. The solid, 
dotted, dashed, long dashed, and dot-dashed lines correspond to Fibonacci, Thue-Morse, pe- 
riod-doubling, Rudin-Shapiro, and random chains, respectively, (a) Initial condition ip n (0) = 
and dijj n (0)/dt = 5 n ,n ; (b) initial condition ip = 5 n ,n and dip n (0)/dt = 0. The chain length 
is N = 10946 for the Fibonacci chain and N = 8192 for the Thue-Morse, period-doubling, 
Rudin-Shapiro, and random chains. 

FIG. 4. The integrated distribution function of frequency (IDFF) as a function of uj at low fre- 
quency for different chains. The solid, dotted, dashed, long dashed, and dot-dashed lines correspond 
to the Fibonacci, Thue-Morse, period-doubling, Rudin-Shapiro chains, and the incommensurate 
FK chain at V < V c , respectively. The results are obtained by directly diagonalizing chains of 
finite length. The chain length is N = 1597 for the Fibonacci and the FK chain, and N = 2048 for 
the Thue-Morse, period-doubling, and Rudin-Shapiro chains. 

FIG. 5. (a) The distribution function of frequency D as a function of u at low frequency for the 
incommensurate FK at V = V c . (b) the time evolution of variance cr 2 (t) in the incommensurate 
FK chain at V = V c for different lengths. The initial condition is ^ n (0) = and dip n (0)/dt = 5 n:Tl0 . 
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